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Abstract 

Two flavor Nambu-Jona-Lasinio model with N components is studied in curved 
space time at finite temperature and density in the leading 1/N expansion. In four 
space time dimension the model exhibits first order phase transition for positive 
curvature. Whereas in flat space an increase in temperature induces a second order 
phase transition even at finite density, in curved space the transition becomes first 
order. We obtain the phase boundary in density-temperature plane and exhibit its 
behavior with changing curvature. A three dimensional plot of the phase boundary 
in /i-T-R space is drawn. 
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1 Introduction 



Spontaneous symmetry breaking (SSB) is one of the most important concepts of 
all unified gauge theories i.e the underlying gauge symmetry is larger than that of 
the vacuum. Of particular interest is the expectation that at high temperatures, 
symmetries that are spontaneously broken today are restored and that during 
the evolution, the universe passed through a series of phase transition evolving 
from a higher symmetric phase to a lower symmetric phase associated with the 
spontaneous breakdown of gauge or perhaps global symmetry. It is very instructive 
to investigate how the phase transition takes place in quantum field theory in 
the environment of the early universe where the temperature, density, primordial 
external electromagnetic field and external gravity may all play important role 

i- 

In the standard grand unified theories, the Higgs field plays a most funda- 
mental role. It is an elementary scalar field with self coupling and has Yukawa 
couplings with all other fundamental fields-fermions and gauge bosons giving them 
masses through its non- vanishing vacuum expectation value. There also exist the- 
ories such as technicolour model of electro-weak theory where the Higgs field 
may be considered as a composite field of some more fundamental fermion fields. 
QCD is an another example of quantum field theory which is invariant under chi- 
ral transformations at the lagrangian level in the absence of quark mass matrix. 
However, the dynamics of the QCD are expected to be such that chiral symmetry 
is dynamically broken with the vacuum state having a nonzero quark-antiquark 
condensate (0 0) and the Goldstone theorem then requries the existence of 
approximately massless pseudo-scalar mesons in the hadron spectrum. To study 
chiral phase transition in QCD we need a non-perturbative treatment and a par- 
ticularly attractive framework to study such systems is the Nambu-Jona -Lasinio 
(NJL) model 0]. Linear sigma model is another such model that is amenable to 
such a study. The NJL model ( for a recent review see [Q and references therein 
) provides a useful framework for studying dynamical chiral symmetry breaking 
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non-perturbatively in the 1/N expansion, where starting from massless quarks, a 
vacuum condensate {qq) 7^ appears and a dynamical fermion mass is generated 
with the breakdown of chiral symmetry in the ground state. Further, NJL model 
can also be considered as a prototype of composite Higgs model and we will use 
it in the present study to investigate the breakdown of higher symmetry in the 
early universe. 

The NJL model four fermion theory in flat space-time in arbitrary dimen- 
sions has been studied using the 1/N expansion method. It is shown that 
chiral sjTumetry is restored in the theory under consideration for sufficiently high 
temperature or chemical potential. It is found that for space-time dimensions 
2 < D < 3 both a first order and second order phase transition occur depending 
on the value of the four fermion coupling while for 3 < D < 4 only the second or- 
der phase transition exists. Recently an investigation of the chiral phase structure 
of the four fermion theory in curved space-time has been made [6-13]. In partic- 
ular the existence of a first order phase transition induced by curvature from the 
chiral symmetric phase to the chiral non-symmetric one has been shown. Chiral 



symmetry structure with nontrivial topoligies |TT], |T3|, |TJ] under the combined ac- 
tion of magnetic fields and curved space-time has also been studied |jl5|. In the 
present paper we will investigate the D=4 NJL model in curved space-time at fi- 
nite temperature, density and curvature using the schwinger proper time method 
[0. We will also take into account the effect of temperature and density on the 
leading contribution on curvature. Linear Sigma model will be topic of our sub- 
sequent study. The existence of first order phase transition at finite temperature 
and curvature is confirmed and the phase diagram in the curvature, temperature 
and density plane is obtained. 
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2 Four Fermion Theory in Curved Space-time 



The Nambu Jona Lasinio model in curved space time is defined by the action 

S= fd^x /(^ [ I V^7^(x) Vm ^ + ^( {^^f + 75 T^ljf) ] (1) 
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where g is the determinant of the space time metric ■,'^^{x) the Dirac matrix in 
curved space-time, \J the covariant derivative of the fermion field ip and N is 
the number of colours, we take the number of flavours to be two. We work in the 
scheme of the 1/N expansion and perform our calculations in the leading order 
of the expansion. For practical purposes it is more convenient to introduce the 
auxiliary fields a and tt and consider the equivalent action. 

S = j d''x^{^)[ii,r{^)y,i'-^{cr^ + ^^)-i^{<y + ^l^rT:)i^] (2) 

Replacing a and vf by the solutions of the Euler-Lagrangian equations arising 
from (2) we reproduce the action (1). This is because the fields a and tt are not 
independent degrees of freedom in (2) and the Euler-Langrangion equations for a 
and vf are infact constraint equations which fix a and tt given ip and ip. 

If a non vanishing vacuum expectation value is assigned to the auxiliary field 
(T, then there appears a mass term for the fermion field ip and the discrete chiral 
symmetry is eventually broken. 

The effective potential (with N factored out) in the leading order of the 1/N 
expansion is then given by 

V{a, = ^ (c^^ + TT^) + z Tr In {x \ {i-f^{x) V^. -(^ + ^ 75 ??) I x) (3) 
In equation (3) the variables a and vr are regarded as constant. Using the 



schwinger proper time method [T^ we rewrite the second term on the right hand 
side of equation (4) as 



V{a, tt) = 77T-(^^ + '^^) + iTrlnS (x, x; s) 
2a 



(4) 

s=a-\- 275 TT 
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where the Green function S defined by 

S{x, y; s) = I {ij^{x) V/. v) 
is the solution of the equation 

{ij^{x) V/. -s)S{x, y; s) = , ^ 6^(x 



y) 



(5) 



(6) 



Thus the effective potential is described by the two point Green's function S(x,x;s) 
of the massive free fermion in curved spacetime. Using the Green function ob- 
tained in the approximation of keeping only linear terms in the curvature, Inagaki, 
Muta and Odintsov H, obtained the effective potential which in D-dimensions 



read as follows: 
y(a,0) 



iTr I as 

2A io 

R 



(2vr) 
1 



Rf,,uk^k%rka + s) 



,2^3 



-YJ'^'^Rcdaiik^Jj^— 



,2^2J 



(7) 



where J"" = 4[7",7 ] and latin indices are vierbein indices. 

The effective potential is divergent in two and four dimensions and is fi- 
nite in three dimensions in the leading 1/N expansion. The four fermion the- 
ory is renormalizable in 2-D fiat space. Using the renormalization condition 



where M is the renormalization scale. The renormalised 



=M 



effective potential in D dimensions is given by 
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We now obtain the four dimensional limit of the NJL model in the MS renor- 
malization scheme given by 

= 2X^M^ -i^(l + 3^-4vr-37)(^) 
1 ^ni^)^---lnAiT + -f) i^j)^ 



^ ' ;n(^)2-l-/n47r + 7) (^)' (9) 



Alternatively one could regularize the divergent part by cutting off the momentum 
integral at finite cutoff A . This gives 

' + + (10) 



(47r)2 6 ^ ' A2 + a2 

The two expression can be shown to be equivalent after carrying out the renormal- 
ization of the coupling constant. The ground state of the theory is determined by 
the minimum of the effective potential (10) namely, by solving the gap equation 



da 



(11) 



For Ar > Ac = 1+3^11-3 7 ' minimum of the effective potential is located at the 
nonvanishing cr, the chiral symmetry is broken down dynamically and a dynamical 
fermion mass is generated. At the critical point the effective potential has two 
degenerate local minima obtained by putting 

y(a, 0) = 1/(m, 0) = (12) 

The solution of the gap equation and the value of the critical curvature can 
be obtained numerically. First, we fix the coupling constant \r greater than the 
critical value Ac corresponding to the broken symmetric phase. To study the phase 
structure in curved space-time we evaluate the effective potential (10) by varying 
the space-time curvature. We see from Fig. la that the chiral symmetry is restored 
as R is increased for a fixed A. The phase transition induced by curvature is of first 
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order as can be seen from Fig. lb. Similar results were obtained by the authors of 
ref. [0 who showed that a first order phase transition is induced by curvature in 
3-dimensional as well. Phase structure of NJL model in curved space time with 
non-trivial topology has been discussed in the literature. It has been shown that 
the combined effect of topology and curvature on phase transition from chirally 
symmetric to chirally broken phase is such as to make the transition first order 
with the growth in curvature as compared to a second order phase transition at 
zero curvature [1^ . 



3 Phase Structure at Finite Temperature and Chemical 
Potential with Curvature 

Four fermion theories in flat space at finite temperature T and chemical poten- 
tial /X in arbitrary dimension has been investigated [5] in the leading order of 
the 1/N expansion. The theory under consideration is renormalisable below four 
dimensions and gives an insight into the phase structure of the theory in four 
dimensions. Following the standard procedure of the Matsubara frequency sums 
[12] we calculate the effective potential in our theory in the leading order of the 
1/N expansion by replacing the phase space integral in eqs (7) by 

'^22 77^ (13) 

n=0 



(27r)4 /3 (27r)3 



where the four momentum k^ is given by 

F = - /i, k) (14) 

and the discrete variable Wn is given by (^"+^^'' 



d^k 



R 



(27r)4^P 



12 (F 



r,2^2 



R 



k^ 



6 (P-s2)3J P^qJ^ (27r)3^P-s 



^ / sds 



d^k 
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R 1 R -. , , 

~12(A;2-s2)2 + 6"(p_s2)3J ^ ^) 

V^{a, 0) = Vo{a) + Vii{a) + V/s{a) (16) 

where 



(i^k , 1 
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Performing integration over s, we get 



2/3 J (27r)3 ^fr'oo " " 24 /? J (27r)3 

~ 1 i R r d^k 

oo 2 

„£o [KTWT^ ^^^^ 

where — k^ + and throwing away the /3 independent terms and carrying out 
the summation over n by the standard techniques we get 

Vp(a) = ^ / k''dk[ln{l + e-/3[V(fc^+-^)-Ml + ^ ^ -/x] 
Ra^ r k^dk 1 

The effective potential ^^((7, 0) at finite temperature, density and curvature is 
then given by 

^''(^^,0) = (l + 3/n47r-37) ^2 

(Zn— - - - Zn47r + 7)(7 



87r2 ' M2 2 



(in — - 1 - /n47r + 7)(j- 
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487r2 ' M2 



-2 



Ra^ r k^dk 1 ^ 

At high densities and low temperatures, the matter is degenerate and the 
fermion distribution function can be approximated by the step function. The 
integral in (20) can be carried out analytically and we obtain 



^ - 2AM^-4^(^ + ^^"^^-^^)m^ 
n [in — /n47r + 7 — - 

[In — - — 1 — /n 47r + 7j 



48M%2 ^ M2 " M2 



37r2 M4 ^ 4 



+-a^ln + 



r2^ 



a " 12M%2 ^ ^ 



+/n (^-^^^^ ^) ] (21) 

But at high temperatures and low densities relevant to the study of phase 
transition in the early universe the integrals in eqs (20) can be done analytically 
by the usual tecniques Ijl^ and we obtain 



— = — — r l + 3 n47r-37 — 



(^'^TT? - 1 - /n47r + 7) 



48M2 7r2 ' M2 " M2 

TT — 1 — ——Ina p H — cj 



7r2/34M4^360 24 32 32 



where c=| + 2 IriA tt — 27. We plot the effective potential(20) as a function of 
a for different values of R at fixed temperature and density. We observe from 
Fig. 2 that at finite temperature symmetry is restored at lower values of curvature 
as compared to the zero temperature case. The solution of the gap equation 
corresponds to the dynamical mass of the fermion and is plotted in Fig. 3a at zero 
chemical potential as a function of curvature for different values of temperature. 
We observe that with the increase in temperature, chiral symmetry is restored at 
lower values of curvature. In Fig.3b, we show that solution of the gap equation as 
a function of chemical potential for different values of curvature at T=0. We again 
find that the symmetry is restored at lower density (lower chemical potential) as 
the curvature increases. 



4 Conclusion 

In this paper we have investigated the phase structure of NJL model in curved 
space time using 1/N expansion and working in the linear curvature approxima- 
tion. But in one discussion we considered rather large values of curvature and it 
would appear that the results obtained may not be reliable. However, it has been 
pointed out 1^ that since terms quadratic and/or higher in R are not divergent, 
Fi? terms are expected to be relatively suppressed compared to linear terms. It 
has also been demonstrated in de-sitter space [0, rT| and Einstein universe |T^] 



that the results obtained in linear approximation are indeed the exact results. 
A distinguishing feature of our investigation is that in the presence of external 
gravity with positive R, the chiral symmetry restoring phase transition is first 
order even as in flat space-time, the transition is second order with temperature 
and becomes first order in the presence of density. In Fig. 4 we have shown the 
phase boundary in temperature, chemical potential plane for different values of 
curvature. We notice that with the increase in R, phase transition takes place at 
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lower temperature and density. In Fig. 5 we have shown the phase boundary curve 
in temperature, chemical potential plane. 
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Figure 1: Behaviour of Effective Potential as a function of for different values of 
the curvature for A > Acr(A = 10). The curves a, b, c, d, e and f are for R = 16, 13, 12, 
8, and -2 respectively. 
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Figure 3: Behaviour of effective potential as a function of cr/M for different values of 
the curvature at fixed temperature T/M =0.5. The curves a, b, c, d, e and f are for 
R/M = 10, 6.5, 4, 2, 0, -1 respectively. 
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Figure 4: Sigma vs curvature for different values of temperature. The curves a, b, c, d, 
e and f are for = 0, 0.2, 0.3, 0.4, 0.5, and 0.6 respectively. 
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Figure 5: Sigma vs chemical potential for different values of curvature. The curves a, 
b, c, d and e are for R/M^ = 0, 2, 4, 6, 8 and 10 respectively. 
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Figure 6: Temperature vs ehemical pot. for different values of eurvature. The curves a, 
b, c, d, e and f are for R/M^ = 0, 2, 4, 6, 8 and 10 respectively. 
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gure 7: Phase diagram in the temperature, curvature and chemical potential plane. 



19 



